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四 阶 拟 线性 椭圆 方程 的 有 限 元 误差 估计 * 
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摘 E 针对 一 类 四 阶 拟 线性 椭圆 方程 ， 本 文 给 出 了 它 的 协调 有 限 元 逼近 。 当 网 格 参数 疡 足够 小 时 ， 得 到 
了 有 限 元 逼近 解 与 真 解 之 间 的 误差 估计 ， 并 且 这 些 误差 估计 是 最 优 的 。 最 后 ， 通 过 数值 实验 验证 
了 理论 分 析 的 准确 性 。 证 明 方法 可 以 类 似 地 应 用 到 某 些 二 阶 拟 线性 椭圆 方程 的 有 限 元 逼近 。 
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四 阶 椭圆 方程 是 一 重要 的 数学 模型 ， 常 见于 板 问题 中 。 因 此 ， 不 论 是 理论 分 析 还 是 数值 分 
析 都 有 很 好 的 应 用 价值 。 从 上 世纪 70 年 代 开 始 ， 有 大 量 的 学 者 都 在 研究 四 阶 问题 ， 特 别 是 在 
有 限 元 的 数值 分 析 方面 。 由 于 Lagrange 元 素 不 能 直接 用 来 求解 四 阶 问题 ， 才 导致 了 非 协 调 有 
限 元 和 混合 有 限 元 的 出 现 上 3。 然 而 ， 上 述 文 献 主要 研究 的 是 标准 的 双 调 和 方程 ， 而 对 于 四 阶 
拟 线性 椭圆 问题 的 研究 并 不 多 ， 据 我 们 所 知 ， 仅 仅 有 Karchevskii 等 的 两 篇 文章 67]。 在 这 两 篇 
文章 中 ， 作 者 分 别 采 用 混合 有 限 元 方法 和 协调 有 限 元 方法 来 处 理 四 阶 拟 线性 问题 。 但 是 对 非 线 
性 项 有 较 强 的 假设 。 他 们 要 求 非 线 性 算 子 满足 强 单调 性 条 件 并 且 是 Lipschitz 连续 的 ， 在 这 种 
情况 下 ， 拟 线性 问题 弱 解 的 存在 唯一 性 可 以 立即 得 到 ， 并 且 真 解 和 逼近 解 的 有 限 元 误差 估计 可 
以 较 容 易 的 得 到 证 明 。 本 文 将 讨论 一 类 四 阶 拟 线性 椭圆 问题 的 协调 有 限 元 逼近 。 在 非 线性 项 有 
较 弱 的 假设 下 ， 我 们 得 到 了 具有 最 优 阶 的 有 限 元 误差 估计 。 最 后 通过 数值 实验 验证 了 理论 分 析 
的 结果 。 

在 本 文 ， 我 们 约定 记号 c 表 示 不 依赖 于 户 的 正常 数 ， 但 是 它 可 能 为 不 同 的 值 。 


2 ”四 阶 拟 线性 椭圆 问题 


考虑 如 下 的 四 阶 拟 线性 椭圆 方程 
A(a(u)Au) =f, 在 9 内 ， 
Bu (1) 
u = Bn = 0, 在 02 E, 


其 中 Q C R2 为 一 有 界 光 滑 的 单 连通 区 域 ，m 为 边界 9Q 的 单位 外 法 线 方向 。 假 定 问题 (1) 的 主 
部 满足 椭圆 性 条 件 ， 即 存在 某 一 正常 数 aal > 0 使 得 对 任意 的 ue€E V = HEN) A 
a(u) > a. (2) 
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我 们 还 假定 a(w) 满足 局 部 Lipschitz 条 件 ， 即 对 任意 的 r > 0， 当 lallv+jalv<r 时 ， 有 
la(w1) — a(ue)| < Llui — uel, (3) 


其 中 工 > 0 依赖 于 +。 给 定 f EV =H), BAAT) 的 变 分 形式 为 : KRuev, 
使 得 
A(u,u,v) =(f,v), vveV, (4) 


其 中 
A(u, vu) = | a(u)dvdwde. 
Q 


对 任意 的 ue V， 定 义 V 中 的 范 数 为 


luv = ( /avpan) ， 


则 根据 (2) 有 A(w,w,w) > ailu]. Kit, WR u eV 满足 问题 (4)， 则 有 
luly < Iglu (5) 


定理 1 在 (2) 和 (3) 的 假设 下 ， 问 题 (4) 至 少 存在 一 个 弱 解 we V， 并 且 满 足 (5)。 
证 明 由 于 V 是 自 反 的 ， 那么 存在 一 组 标准 正 交 基 {®;}%| EV PARARE. W Vm Ho 
Bl Om 所 张 成 的 空间 ， 定 义 问 题 (4) 的 Galerkin EAH un 满足 


A(Um, Um, Bi) = (f, Bi), i=1,---,m, (6) 
JEU um 的 存在 性 可 由 Brouwer 不 动 点 定理 得 到 。 根 据 (5)， 有 
lunlv < 二 lw: 


因此 ， 存 在 {um} S 的 子 序 列 RAWEA {um} L) Mu Ee V， 使 得 当 m 一 +co 时 ， 
Ai um HEV PRM Fu, ELO PRR u, FAQ PIL FOAM Fu. FRR 
ATWEBA u HAE [al Bl (4) 的 弱 解 。 


|A(Um, Um, ®;) =. A(u, u, ®;)| 


IA 


|A(um, Um, Bi) 一 A(u, Um, Bi)| + |A(u, Um, Bi) 一 A(u, u, ;)| 
= h+. (7) 


ER, RNA 
h= j alu)A(um — u)A®idz > 0, 4 mm 一 oo. 
Q 


根据 (3)， 我 们 还 有 
i= J (alum) — a(u))Au,A®;dx 
Q 
< Lum — ullrllum|lv|Pilly +0, 4 m— +00. 


因此 ， 在 问题 (6) 中 取 极 限 并 注意 到 ®; 的 性 质 ， 我 们 就 证 明了 v 是 问题 (4) 的 一 个 弱 解 。 
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3 “有 限 元 逼近 
WT, ARENA Ae (或 四 边 形 ) 剖 分 ， 其 中 参数 0 < h < 1。 定 义 有 限 元 空间 


Jie C1(Q)N H?(Q), ulre Ps, VT ETh} 三 角形 前 分 ， 
{u € C1(Q)N H7(Q), u |E Q3, YT € Th} WAER. 
WV, = VAWr MV AV WA BCT lal. WAHE (4) 的 有 限 元 逼近 形式 为 : Kun E 


Va» 18 
A(up, Un; Vh) = (f, vh), V un E Vh. (8) 


定义 算 子 工 è Va -> Va» 满足 
Alun, Tun, vn) = (f, un), V Un € Vn. (9) 
WES 
D = {un € Vr, llually < Hfl} 
LA -一 Ql 


APERT AA DED 的 一 个 连续 映射 ， 则 由 Brouwer 不 动 点 定理 有 : 
定理 2 在 (2) 的 假设 下 ， 离 散 问题 (8) 至 少 存在 一 个 解 u © Vie 


4 误差 估计 
为 了 得 到 有 限 元 误差 估计 ， 我 们 要 求 a(w) 关于 久 是 二 次 连续 可 微 的， 并 且 满 足 局 部 性 条 件 
a(u)<a2, Vue B(R), (10) 


其 中 B(R) = {v € V, llully < R}, az > 0 依赖 于 R。 下 面 我 们 再 引入 一 些 有 限 元 的 基本 假 
wt, 
(Aj): 


A = < Airs 
inf, lu valla < ch~ lulas 


HHO<s <2. 
(A2): 
llenlly < ch~*llvallo, Wun E Vh, 


SEF |lunllo = lwallzz。 
定理 3 在 (2),(3),(10) 和 (41)-(42) 的 假设 下 ， 如 果 w e HA4NV 和 wh E Vn RA i 
题 (4) 和 问题 (8) 的 解 ， 那 么 对 足够 小 的 h， 我 们 有 如 下 的 最 优 阶 误差 估计 


ju ~ ually < ch?,  |lu—unllo < chê. 
证 明 在 问题 (4) PR = vs， 并 与 问题 (8) 相 减 得 
A(u, u, Un) 一 A(un, Uh, va) =0, Vun € Va. (11) 


则 
A(up,U— Un, Un) = Alun, U, Un) ~ Alu, u, Un). 
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因此 ， 根 据 (2), (10) 和 (11) 有 
ollu — all? < Alur, u — up, u— un) 
= A(un,U— Uh, U — Up) + A(Uh, U — Un, Uh — Un) 


< a2|\u — ually lu — vally + Alur, u, vn — un) — Alu, u, Uh — Un) 


az|lu — ually lju — vally + T. (12) 


下 面 估计 六， 
1 = f (elu) ~a(w)dud(oy -ude < L Ju — us duid(on — wladz 
和 Lull, lju — uallollua — ally < ch ?lv — uallo(Ilu — wnllo + lu — vallo) 
< ch ?lu — wallg + ch? lu — uallollull za. 
将 上 式 代 入 (12) 得 
orlu — wally < allu — uallvilu — vally + ch? |lu — walle + ch ?le — ual 
< ch? l|u— ually + eh- ?lu — unl. 


应 用 Young TER, A 


lu — wnlly < cht + ch ?lu — uals. (13) 


下 面 估计 |u- wnllo。 考 虑 问题 (1) 的 导 算 子 方程 


Lw =p, ENA, 
(14) 
w= = =0, £00 上 ， 
其 中 
Lw = A(a(u)Aw + Dya(u)Auw). 
问题 (14) 的 对 偶 方程 为 
L*w = q, 在 Q A, 
(15) 
we ae =0, #ank, 
其 中 


L*w = A(a(u)Aw) + Dya(u)AuAw. 


对 任意 的 ge L?(Q), AE (15) 总 存在 一 个 解 w e HANV, FH lwla < cllqllo。 又 由 有 限 元 
的 插值 理论 ， 存 在 wn e Vi 使 得 


llw — wally < ch?l|wl|zr4 < ch?llallo- 
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在 问题 (15) P, Sqau-—u, WA 


lu url = (L*w,u— ur) = (L(u — ua), w) 


| 


TEONE — un)Aw + Dya(u)Au(u — up)Awdz 
= A(u,u,w) — A(u, un, w) +f Dya(u)Au(u — un)Awdz 
Q 
= A(u,u,w) 一 Alun, uh, w) + f (alun) — a(u)) Aur Awdz 
2 

+ [ Duawau(u — un) Awdz 

= A(u,u,w — wp) — Alun, Un, wW — wh) + f D.,a(ū)AunrAw(un — u)dx 
Q 


+f D,za(u)AuAw(u — up)dz 
2 


IA 


Alur, Un — u, w — wp )| + |A(u, u, w — wr) — Alun, u, w — wh )| 


+ f D,a(ū)(Aun — Au)Aw(un — u)dz 


+ f (Duat) — Dya(u))AuAw(u, — ujde 
Q 


< ellu ~ urllv lw — wallv + ellu — urlizs llu — unliv lwll 
< ch?|lu — ually llu — wallo + cllu — unllze llu — ually ilu — uallo- 
则 由 Young 不 等 式 有 
lu — urlo < ch? |lu — ually + cllu —uallzsilu — uallv 


ch? ||u — ually + cllu — wally llu — walle 


lA 


IA 


ch? ||u — ually + cllu — ually + ellu — uallo. 
取 s = 二， 我 们 得 到 
llu — unllo < ep2llu — ually + cllu — ually. 
将 (16) (RA (13) 得 
lju — wnlly < chf + ch? lu — ually + ch? lu ~ ually. 
因此 ， 对 足够 小 的 h”， 有 
lu — unlly < cht + ch ?lu — ually- 


# llu — unlly < ch， 则 
lu — unlly < chê + ch?|lu — ually- 


那么 对 足够 小 的 有 有 |u — ually < ch?。 将 上 式 代入 (16) 即 得 ||w — unllo < ch4。 
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(16) 
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5 ”数值 实验 
Balu) = 1+ u? WA (2), (3) 和 (10)， 因 此 ， 我 们 考虑 


A((l+u7)Au) =f, ENA, 
=z 
其 中 Q = [0,1] x [0,1] WAH u = sin? rsin? y, MARRU u WER KAUR RE. 

我 们 考虑 四 边 形 前 分 ， 有 限 元 素 为 Bogner-Fox-Schmidt 7t (Q3 元 )。 假 设 z 方 向 与 y 方 向 的 
剖 分 数 相同 ， 并 都 记 为 及。 表 1 分 别 给 出 了 在 不 同 的 网 格 下 ， 有 限 元 求解 的 HH? 误差 和 52 误 
差 。 我 们 可 以 看 出 ， 当 网 格 不 断 加 密 时 ，H? 误差 和 L? 误差 越 来 越 小 ， 并 且 误 差 阶 数 也 与 我 们 
的 理论 分 析 较 吻合 。 


0, 在 OQ E, 


u 


Al: 不 同 网 格 下 的 误差 及 其 阶 数 


0.282592 


8.154853 x 1073 
4.352941 x 1074 
3.241255 x 1075 


2.102684 x 1078 


4.975746 x 10-3 
| 4.00 | 1.244261 x 1073 
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Finite Element Error Estimate of Fourth Order 
Quasilinear Elliptic Equations 
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Abstract: For a class of fourth order quasilinear elliptic equations, the conforming finite element 
approximation is presented. If the parameter h of the grid is sufficiently small, the error estimates 
between the finite element approximation solution and the true solution are obtained. Moreover, these 
estimates are optimal. Finally, the numerical experiment is given to verify the correctness of the 
theoretical analysis. The proof method can be extended to the finite element approximation of some 
second order elliptic equations. 

Keywords: fourth order quasilinear elliptic equation; finite element approximation; error estimate 
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